Controlling ultracold Rydberg atoms in the quantum regime 
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Abstract 

We discuss the properties of Rydberg atoms in a magnetic Ioffe-Pritchard trap being commonly 
used in ultracold atomic physics experiments. The Hamiltonian is derived and it is demonstrated 
how tight traps alter the coupling of the atom to the magnetic field. We solve the underlying 
Schrodinger equation of the system within a given n-manifold and show that for a sufficiently 
large Ioffe field strength the 2n 2 -dimensional system of coupled Schrodinger equations decays into 
several decoupled multicomponent equations governing the center of mass motion. An analysis of 
the fully quantized center of mass and electronic states is undertaken. In particular, we discuss the 
situation of tight center of mass confinement outlining the procedure to generate a low-dimensional 
ultracold Rydberg gas. 

PACS numbers: 32.60.+i,33.55.Be,32.10.Dk,33.80.Ps 



1 



Rydberg atoms possess remarkable properties. Although being electronically highly ex- 
cited they can possess lifetimes of the order of milliseconds. Moreover, due to the large 
displacement of the valence electron and the atomic core they are highly susceptible to elec- 
tric fields and therefore easily polarizable. The latter is impressively shown in ultracold 
Rydberg gases [lj where the mutual multipole interaction of two Rydberg atoms leads to 
a number of intriguing many body effects, such as resonant Forster transitions [2] and the 
dipole blockade mechanism 13, 14|. Their tunable electric dipole moment makes them also 

n 

interesting candidates for the realization of a two-qubit quantum gate |5J which is a crucial 
ingredient for the realization of quantum information processing schemes. However, a pre- 
requisite for realizing the latter in a coherent fashion is the availability of a well-controlled 
environment in which the Rydberg- Rydberg interaction takes place: Firstly one has to assure 
the Rydberg atoms to be individually addressable. Secondly, since the mutual interaction 
strength strongly depends on their relative displacement, a suitable way to spatially arrange 
the atoms has to be found. 

There exist several proposals for building traps for Rydberg atoms suggesting the use 
of electric jfj, optical 01 and magnetic fields [Sj . The experimental feasibility of magnetic 
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trapping has been shown recently by employing very strong magnetic fields [9j|. Trapping in 
the quantum regime, which is indispensable to gain precise control over the atomic motion, 
however, could not yet be demonstrated. Properties of such quantized Rydberg states that 
emerge in a magnetic quadrupole trap have been theoretically discussed in Refs. [§]. The 
authors report on a strong confinement of the atomic center of mass (cm.) with trapping 
frequencies of the order of 10 kHz. At the same time the electronic structure is, to a large 
extent, unchanged compared to the field-free case. Although the 3D magnetic quadrupole 
field is omnipresent in ultracold atom experiments it has one drawback. Due to the point of 
zero field strength at its center Majorana transitions are likely to happen thereby leading to 
loss from the trap. In Ref. js| it was shown that by increasing the total angular momentum 
of the atom this problem can be overcome. However, for practical purposes it is desirable 
to have trapping at small or even zero center of mass angular momentum since this is the 
regime in which trapped ground state atoms are usually prepared. 

In this letter we show that trapped and controllable Rydberg states can be achieved in a 
Ioffe-Pritchard (IP) trap. We discuss how the large size of the Rydberg atom modifies the 
coupling to the magnetic field in comparison to ground state atoms and demonstrate the 



2 



feasibility of generating low dimensional Rydberg gases. Using the IP configuration is of 
particular relevance since it is usually employed for preparing ultracold atomic gases which 
can serve as the initial state for the production of trapped ultracold Rydberg atoms. 

In the two-body picture the mutual interaction of the valence electron (particle 1) and 
the remaining core (particle 2) of the Rydberg atom can be modeled by an effective potential 
that exhibits short range properties describing the electron-core scattering and a long range 
Coulombic behaviour. States with large electronic angular momenta I, which represent 
the focus of the present investigation, probe almost exclusively the Coulombic tail of this 
potential. We do not account for relativistic effects such as spin-orbit coupling which for 
large n, I are negligibly small compared to the energy shift due to the magnetic field of the 
IP trap (sl Q| • The interaction of the electronic and core charge with the magnetic field are 
taken into account via the minimal coupling. Following the above arguments and including 
the coupling of the nuclear and electronic magnetic moments to the external field our initial 
Hamiltonian reads 

ff,ni. = ^[pi-eA( ri )] ! + 5 -L[p 2 + eA(r 2 )r 

+V(\r, - r 2 |) - Ml ■ B(r,) - M2 ■ B(r 2 ). (1) 

The IP field configuration is given by B(r) = Be z + Bi in (r) + B qua d(r) with the 
linear component Bi in (r) = G [xe x — ye y ] and the quadratic component B quad (r) = 
Q [2z(xe x + ye y ) + (x 2 + y 2 — 2z 2 )e z )}. In the following we consider the case of a domi- 
nant linear component and neglect B qua d(r). We therefore encounter a two-dimensional 
quadrupole field which is translationally invariant along the z— axis. Our corresponding vec- 
tor potential reads A(r) = A c (r) + Ai in (r) with A c (r) = -| [xe y — ye x ] and A lin (r) = Gxye z . 
Next let us insert these expressions into the Hamiltonian and introduce relative r and 
cm. coordinates R with their respective momenta p and P. In the absence of the external 
field one arrives at a decoupled cm. and electronic motion. However, already in the presence 



of a homogeneous field this is not the case 



ll| and consequently terms coupling the cm. and 



internal motion emerge. To simplify the latter in case of our inhomogeneous magnetic field 
configuration we apply the unitary transformation U = exp e z x r ■ R] which eliminates 
couplings of the cm. and the relative motion generated by the Ioffe field (atomic units are 
used throughout except when stated otherwise) Furthermore we neglect the diamag- 
netic terms, which is a very good approximation for laboratory field strengths and not too 



high principal quantum numbers n (see ref. |8|), and keep only the leading order terms with 
respect to the inverse masses 

p2 

H 1P = H A + — -n 2 - B(R) - fii ■ B(R + r) + A lin (R + r) ■ p. (2) 

2IVI2 

Here Ha — \ — - is the Hamiltonian of a hydrogen atom possessing the eigenfunctions 
\n, I, mi, m s ) and energies E n = —^n~ 2 . The following two terms of Hip describe the cm. mo- 
tion of a point particle possessing the magnetic moment u in the presence of the field B. 
This system has been thoroughly investigated in Refs. Since the magnetic moments 

are connected to the corresponding spins S and £ according to fi 1 = — S and fj, 2 = — ^jH, 
with being the nuclear g-factor, we neglect the term involving [i 2 in the following. The 
last two terms of Hip couple the electronic and cm. dynamics mediated by a spin-field 
and motionally induced coupling. We remark that the Hamiltonian (0) commutes with the 
z-component of the linear cm. momentum P z being a direct consequence of the above- 
mentioned translational invariance of the system along the z-axis. Hence the longitudinal 
motion can be integrated out employing plane waves \K Z ) = exp (—iK z Z). In order to solve 
the eigenvalue problem of the resulting Hamiltonian, that depends on five spatial degrees of 
freedom, we assume the magnetic field not to couple adjacent n- manifolds. Estimating the 
energy level shift caused by the magnetic field according to E zee « Bn this requirement is 
fulfilled M\E n - E n ^\/E zee « (5n 4 ) > 1. In this case we can consider each n-manifold 
separately and represent the Hamiltonian (J2J in the space of the 2n 2 states spanning the 
given n-manifold. Neglecting the constant energy offset E n and introducing scaled cm. co- 
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ordinates (R — > 7^3 R with 7 = GM 2 ) while scaling the energy unit via e sca i e = 73 J Mi we 
eventually arrive at the working Hamiltonian 

P 2 + P 2 

n = y + ^ _ G(X, Y) + 7 -3M 2 H r (3) 

with the effective magnetic field 

G(X, Y) = Xe x - Ye y + ^M^e, (4) 

where /x and 7i r are the 2n 2 -dimensional matrix representations of the operators 1/2 [L + 2S] 
and H r = Ai in (r) • p + Bi in (r) • S, respectively. Here we have introduced the orbital angular 
momentum operator L = r x p. 

The Hamiltonian (JSJ) can be interpreted as follows: The first two terms describe the 
dynamics of a neutral atom in a IP trap under the assumption that the coupling of the atom 



to the field is given by the Zeeman energy E z = 1/2 [L + 2S] • G(X,Y). One encounters 



a similar coupling term also for 'point-like' ground state atoms [13( where, however, the 
generically strong hyperfine coupling leads to the fact that \x is proportional to the total 
spin F. The distinct feature of the Hamiltonian © is the appearance of the last term 
which accounts for the finite size of the Rydberg atom. This term scales according to 
7~3GM2n 2 = 7^n 2 ~ 7^ (r). Since 7"^ can be regarded as a length unit for the cm. wave 
function we find this term to be particularly important if (r) w (R), i.e. if the size of the 
atom and the size of the cm. state become comparable. In a typical macroscopic IP trap l| 
the cm. wave functions are very extended and the energy spacing between the cm. states 
is small. Hence, in this 'classical' regime 7i r has little effect and can be neglected. The 
situation changes for tighter IP traps which can, for example, be realized by so-called atom 
chips Here the micro structured wires allow for the generation of traps for which the 
size of the cm. ground states are of the order of several 100 nm and become therefore 
comparable to the typical size of Rydberg atoms. 

In order to solve the Schrodinger equation belonging to the Hamiltonian (J3J) we employ 
an adiabatic separation of the electronic and the cm. motion. To this end a unitary trans- 
formation U (X, Y) which diagonalizes the last two (matrix) terms of the Hamiltonian, i.e. 
W(X, Y)(n- G(X, Y)+ 1 -lM 2 H r )U(X, Y) = E a (X, Y), is applied. Since U(X, Y) depends 
on the cm. coordinates the transformed kinetic term involves non-adiabatic (off-diagonal) 
coupling terms which we will neglect in the following. We are thereby lead to a set of 2n 2 
decoupled differential equations governing the adiabatic cm. motion within the individual 
two-dimensional energy surfaces E a , i.e. the surfaces E a serve as potentials for the cm. of 
the atom. Figure ^ shows an intersection through a subset of these surfaces for the case of 
87 Rb in a IP trap with a gradient G = 20 T/m and a Ioffe field strength B = 0.01 G. One 
immediately notices the large number of avoided crossings between the adiabatic potential 
surfaces. Here non-adiabatic transitions mediated by the (neglected) off-diagonal coupling 
terms of the kinetic energy are likely to occur. The uppermost surface, however, does not 
exhibit such avoided crossings and is a possible candidate in order to achieve stable trap- 
ping. According to our findings this surface is predominantly formed by the electronic state 
possessing the largest possible orbital angular momentum in the n-manifold under consid- 
eration. The corresponding quantum defects and relativistic corrections are therefore tiny 
which a posteriori justifies their neglect. 
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FIG. 1: Cut through the adiabatic potential surfaces of the n = 30-manifold ( 87 Rb, G = 20 T/m, 
B = 0.01 G). The high density of states is clearly visible. The uppermost potential surface is 
predominantly formed by the atomic state with largest angular momentum. This particular surface 
is clearly separated from the next lower ones whereas all other surfaces exhibit a number of avoided 
crossings (see the magnified view in the inset) at which non-adiabatic inter-surface transitions are 
likely to occur. 

The actual appearance of the potential surfaces depends on both, the gradient and the 
strength of the homogeneous Ioffe field. In figure ^ the latter is comparatively small such 
that for small displacements from the coordinate center its field strength is easily sur- 
passed by that of the gradient field. The different symmetry properties of the both fields 
lead to a rearrangement of the energy surfaces finding its expression in the large num- 
ber of avoided crossings. We now turn to the case where the Ioffe field is large com- 
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pared to that of the gradient field, i.e. r )~^M.2B ^> 1. This however can hold true 
only in the vicinity of the trap center. We assume this region to be sufficiently large 
such that at least a few low-lying cm. states are localized here. In this case the term 
fx - G(X, Y) will dominate the Hamiltonian (j^J. We now diagonalize this term by employing 
the transformation Ud(X, Y) = exp [i(L z + S z )<fi) exp [i(L y + S y )j3] with </> = arctan [Y/X], 
cos/5 = -f-lM 2 B\G(X,Y)\- 1 and sin/3 = — (X 2 + Y 2 )|G(X, Y)| -1 . This yields the adiabatic 
energy surfaces 

E a = Ul(X,Y)n.G(X,Y)U D (X,Y) = ±(L Z +2S S )\G(X,Y)\, (5) 

being characterized by the quantum numbers of L z and S z which are mi and m s , respectively. 
The energetically highest surface is assumed for mi = n — 1 and m s = 1/2. The next lower 
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FIG. 2: Surface plot of the the seven uppermost energy surfaces of the n = 30-manifold ( 87 Rb, 
G = 20T/m B = 0.1 G). Clearly the grouping into three multiplets whose mutual distance is given 

_ 2 

by 7 3M%B is visible. While the uppermost manifold consists only of a single surface the next 
lower ones show an approximate twofold and fourfold degeneracy. A magnified view of them is 
provided in the insets. 

one is twofold degenerate and the following one shows a fourfold degeneracy (see figure |2J). 

2 

The energy gap between these degenerate multiplets is given by AE rs ^~zM 2 B and can 
hence be continuously varied by tuning the Ioffe field strength. In the present regime the 

2 

term •y~3M 2 7i r can be considered as a perturbation since its energetic contribution is much 
smaller than AE. The correction to the uppermost surface is zero whereas 7i r couples the 
surfaces of the energetically lower lying degenerate multiplets. To study the dynamics of 
the multi-component cm. wave function within these coupled potential surfaces constitutes 
a very interesting problem. In the present investigation, however, we will focus exclusively 
on the uppermost non-degenerate surface. 

The explicit knowledge of Ud(X, Y) allows for an analytical calculation of the non- 
adiabatic couplings between any of the potential surfaces arising from the kinetic energy 
term. Our findings show them to be proportional to ^{M 2 B)^ 1 = (AE)^ 1 . For a suf- 
ficiently large Ioffe field strength we can thus safely employ the adiabatic approximation, 
i.e. neglect the non-adiabatic coupling between the uppermost and the next lower sur- 
face. In order to obtain the quantized cm. states we solve the scalar Schrodinger equation 
[l/2(Pj + Py) + E a (X, Y) — €v\ $? U (X, Y) = in the uppermost potential surface which we 
denote by a = 2n 2 . For low cm. excitations the potential surfaces are approximately har- 
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FIG. 3: Expectation value of the radii p 



VX 2 + Y 2 and r of the cm. and electronic wave 



function in the IP trap ( 87 Rb, G = 100 T/m, B = 0.1 G). v labels the cm. quantum states within 
the uppermost adiabatic potential surface. While {p) is increasing (r) remains approximately 
constant at its field free value. For small degrees of cm. excitation the cm. state is even stronger 
localized than the valence electron ((p) < (r)). Moreover, since the electron is found in a high 
angular momentum state its radial uncertainty Ar is small. Thus a scenario where the cm. and 
the electronic wave function do not overlap is possible as sketched in the inset. 

monic and thus the energies are in reasonably good agreement to those of a two-dimensional 
isotropic harmonic oscillator with a n-dependent frequency uj = G 2 BMi ( m a ^ om ^ c units). 
Hence, by choosing high gradients and an appropriate tuning of the Ioffe field strength B 
very tightly confining traps for highly excited atoms can be obtained. Such a situation is 
depicted in figure El where we show the expectation values of the radii p = \/X 2 + Y 2 and r 
for the cm. and electronic wave function versus the degree of excitation of the cm. motion 
v. In the presented case the confinement gives rise to a trap frequency of approximately 
1.4 MHz. In this regime the size of the cm. state characterized by (p) is even smaller than 
the electronic cloud, i.e. the cm. wave function is stronger localized than the valence elec- 
tron. On the other hand the expectation value (r) for the electron remains nearly constant 
possessing the corresponding field free value as the degree of excitation of the cm. increases. 
This indicates that in spite of the strong localization of the cm. the electronic structure 
of the atom is barely changed compared to the field free case. This observation has been 
backed up by calculating further electronic properties, such as the expectation values of L 2 
which also barely differ from their corresponding field free values. As previously indicated 
we find the electron in the highest angular momentum state (I — n — 1) which possesses the 
smallest radial uncertainty Ar for given n. Due to this fact it is possible that the cm. and 
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the electronic wave function may not even overlap (see inset of figure EJ). This novel regime 
opens up the possibility to control Rydberg atoms in the quantum regime and might pave 
the way to study many-body effects in low- dimensional ultracold Rydberg gases [14j|. 

In order to study the latter, ultracold atoms confined in a tight atom chip trap [l5| can 
be transferred into high-/ Rydberg states by imposing suitable optical and radio frequency 
fields (see refs. jlfl 13 an d re fs- therein). Since the electronic structure is barely affected 
even for tight cm. confinement the Rydberg atoms keep their well-known properties such 
as long radiative lifetimes and electric dipole moments. 
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